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Abstract. We obtained N- and Q-band observations of the Apollo-type asteroid 25143 Itokawa during its close 
Earth approach in July 2004 with TIMMI2 at the ESO 3.6 m telescope. Our photometric measurement, in combi- 
^0 ' nation with already published data, allowed us to derive a radiometric effective diameter of 0.32±0.03 km and an 

cal and flattened ellipsoid of the approximate size of 520(±50) x270(±30) x230(±20) m, based on the Kaasalainen 
ly-^ ' et al. 120051 1 shape model. Our studies show that the thermal observations lead to size estimates which are about 

15% smaller than the radar results (Ostro et al. 120051 . slightly outside the stated radar uncertainties of ±10%. We 



Oh 



determined a rather high thermal inertia of 750 J m -2 s~ u ' b K _i . This is an indication for a bare rock dominated 
surface, a thick dust regolith can be excluded as well as a metallic surface. From our data we constructed a 
10.0 /im thermal lightcurve which is nicely matched in amplitude and phase by the shape and spin vector solution 
in combination with our TPM description. The assumed S-type bulk density in combination with radiometric size 
, lead to a total mass estimate of 4.5l 2 'g ■ 10 10 kg. 
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1. Introduction Goldstone, resulting in size, shape, radar albedo and sur- 

& ' face roughness estimates. Kaasalainen et al. i|2003l 12005) 

The Near-Earth asteroid (NEA) 25143 Itokawa determined through lightcurve inversion techniques a high 

(1998 SF36) is the target for the Japanese Hayabusa qua i ity shape and pole solution together with a period of 

(MUSES-C) sample return mission. The spacecraft will p M = 12.13237 ± 0.00008h. Sekiguchi et al. (2003) ob- 

arrive at the asteroid in summer 2005 and hover close served Itokawa at thermal infrared wavelengths and used 

to the surface for about 3 months before it will collect the powerful radiometric technique (e.g. Harris & Lagerros 

surface samples which will be brought back to Earth in to determ i ne a size of 0.35±0.03km and an albedo 

June 2007. of .23(+0.07, -0.05). Miiller et al. (2004) reported on 

Several ground-based observing campaigns took place thermal property studies of Itokawa based on multi-epoch 

during the last years to d erive various properties of thermal infrared photometric data. But their data set did 

Itokawa. E.g., Binzel et al. l|2001f> concluded from vis- no t a ll ow to find a unique solution for the size, thermal 

ible and near-infrared spectroscopic measurements that inertia and roughness. Based on the radar size, they de- 

the spectral characteristics (S(IV)-type) matc h the LL or- r i ve d a thermal inertia between 5 and 10 times that of the 

dinary chondrite class meteorites. Ostro et al. (200H l2TiU5)) Moon, 
report on delay-Doppler images obtained at Arecibo and 

Here, we used all available thermal infrared data to- 



Send offprint requests to: T. G. Miiller gether with own observations, taken during a close ap- 

* Based on observations collected at the European Southern proach in July 2004 (Section EJ). The modelling and the 
Observatory, Chile; ESO, No. 73.C-0772 derivation of thermophysical properties of Itokawa was 
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Table 3. Summary of TIMMI2 observational results for 
asteroid 25143 Itokawa. Note: The Sekiguchi et al. p^B) 
and Delbo |2004) fluxes and errors have been recalculated 
for the true central filter wavelengths, using the corre- 
sponding stellar fluxes. Measurements from April 9, 2001 
(marked with *) were taken under less favorable condi- 
tions in comparison with data from April 8 (M. Delbo, 
priv. comm.). Our new observations are listed under #12 
until #20. 



then performed through the well-established, tested and 
frequently used thermophysical model (TPM) by Lagerros 
((MR)! ITM71 IT3951 see Sects. Hand 0). The results are 
discussed in the context of the already known properties 
of Itokawa (Sect.EJ) and conclusions are drawn in Sect. [7J 

2. Observations and Data Reduction 

We combined the observations from Sekiguchi et al. (P003) 
(data set #1 in Tabled, with the data set from Delbo 
(2004J (data set #2), and our own observations (data set 
#3). All measurements were taken with the TIMMI2 in- 
strument (Kaufl et al. OTMjl at the ESO La Silla 3.6 m 
telescope. Table ^ summarises the observing geometries 
for all 20 measurements. 

A standard chopping and nodding technique was uti- 
lized for all observations to reduce the atmospheric and 
telescope background emission. Chop and nod throws were 
10", respectively. For the imaging observations, a pixel 
scale of 0.2" was chosen, and on source integration times 
were 22min (obs #1), 10-13 min (obs #2 to #11), 7.5 min 
(obs #12 to #19), and 25 min (obs #20). 

Both already published data sets, from Sekiguchi et 
al. H2003f) and from Delbo f2003)l. were re-calibrated using 
the official central filter wavelengths in combination with 
the flux densities of the corresponding stellar models (see 
Tbl.©. 

A more detailed data reduction and calibration was 
performed on our own data set. Unfortunately, the tele- 
scope tracking on 25143 Itokawa was not perfect in the 
2004 observing run. This had the consequence that the 
pipeline-reduced images showed elongated sources (see left 
side of Fig. QJ. Therefore, we took the raw images and 
processed them in a pipeline-like manner together with a 
fixed x-y-shift rate of typically 0.5 to 1.0 pixels/minute. 
A centroid determination in combination with a standard 
shift-and-add technique was not possible due to the faint 
sources which are not always visible on individual raw 
frames. The results of this process are illustrated on the 
right side of Fig. ^ 

We applied standard aperture photometry on the 
tracking corrected images with same apertures on the 
stars and asteroids. The aperture radii were chosen 
with respect to the growth curves (for details see 
e.g. Delbo [2004). All four on-array signatures of the 
sources were used for the flux calibration. Colour 
differences between stars and 25143 Itokawa were 
negligible (about 1-3 %) for the used filters in combi- 
nation with the atmospheric transmission at La Silla and emitted energy, and also the heat conduction into 
(http://www.ls.eso.org/lasilla/sciops/3p6/tiimi/httfi^ and emissivity control 

The observational results are summarised in Tbl. yl the energy balance and thereby the surface temperature. 

The error values of observations #12 to #20 include The thermal inertia in combination with the rotation pe- 
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error calculations of observations #1 to #11 are described 
in the corresponding references. 

3. Thermophysical model description 

We applied the TPM by Lagerros <P<)6i 110071 IT^Sj) to all 
20 measurements from Tbl. |21 to investigate the physical 
and thermal properties of Itokawa. On the large scale, the 
TPM considers the asteroid size, the global shape and spin 
vector and the actual observing and illumination geome- 
try at the time of an observation. On the small micrometer 
scale, the TPM takes into account the reflected, absorbed 



the uncertainties of the calibration star models (3 % in N- 
band and 4% in Q-band), the aperture photometry error 
(Nl: 3%, N2: 2%, N12.9: 2%, Ql: 10-20%) and an error 
for the flat-field residuals (about 3-6 % in N and about 10- 
15% for Q-band, depending on the relative placement of 
the calibrator and the asteroid on the TIMMI2 array). The 



riod and the orientation of the spin vector influence the 
diurnal temperature variations. As a result, the thermal 
inertia is strongly connected to the interpretation of mid- 
IR observations, namely when comparing before and after 
opposition observations at large phase angles with very 
different temperatures of the terminator. Moreover, the 
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Table 1. Summary of TIMMI2 observations of asteroid 25143 Itokawa. The phase angles are negative before opposition 
and positive after. We added data from Sekiguchi et al. I|2003fl and Delbo H2004|l . 
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Table 2. Monochromatic flux densities in [Jy] of the stellar calibrators at the TIMMI2 central filter wavelengths. Note: 
The key wavelengths were taken from jhttp: //www .Is .eso.o rg/lasilla/sciops/3p6/timmi/Filte rsl The model 
fluxes were taken from http : //www. iso . vilspa. esa. es/users/explJ.ib/ISD/wwwcal/ isoprep/cohen/templates/ 
and interpolated to the central filter wavelengths. 
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17.41 
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thermal inertia determines the amplitude of the thermal 
lightcurve for a given aspect angle. The TPM beaming 
model, described by p, the r.m.s. of the surface slopes and 
/, the fraction of the surface covered by craters, accounts 
for the non-isotropic heat radiation, noticeable at phase 
angles close to opposition. But it also influences the shape 
of the spectral energy distribution in the mid-IR. Detailed 
considerations of the T, p and / influences for various ob- 
serving geometries and wavelengths are discussed in e.g., 
Miiller lj2T)Tl2")l or Dotto et al. pTO| . 

An iJy-value of 19.9 (Kaasalainen et al. l2UU3l M. Abe, 
priv. comm.; Sekiguchi et al. I2003f) and a G- value of 0.21 
(Abe et al. I2()()2al I2002b(l was used to describe the visual 
brightness of Itokawa. We assumed a constant emissivity 
of 0.9 at all wavelength. 

The shape- model used here (Kaasalainen et al. [5005 ) 
is an update and refinement of the model presented 



in Kaasalainen et al. f^S3- The model accommodates 
new photometric observations from December 2003 to 
September 2004, as well as some 2001 data additional to 
the 2000-2001 apparition dataset presented in Kaasalainen 
et al. H2003J) . The long time-line of the updated dataset 
allowed accurate period determination for Itokawa, and a 
refined pole and shape estimate. All parts of Itokawa's sur- 
face were well visible during the two apparitions; however, 
the long period precluded fully covered rotational phases 
for single lightcurves. Calibrated photometry allowed the 
determination of Itokawa's solar phase curve for a wide 
range of solar phase angles. The refined rotation parame- 
ters are (3 = -89° ± 5°, A = 330° for the ecliptic latitude 
and longitude of the pole, and P = 12.13237 ± 0.00008 h 
for the sidereal period. Fig. shows equatorial edge-on 
and pole-on images of the shape model. The model agrees 
well with the radar-based one (Ostro et al. 1271775) . 
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Fig. 1. Pipeline-processing in comparison with re-reduced asteroid images. Top: Itokawa measurement #12 (Tbl. 
Nl-filter, 7.5 min integration time, left: pipeline-product, right: the 60 raw images were shifted by 1.0 pixels per minute 
in x-direction and by -1.0 pixels per minute in y-direction and then co-added. Bottom: Itokawa measurement #20 
(Tbl.^), Ql-filter, 25.8 min integration time, left: pipeline-product, right: the 234 raw images were shifted by 0.8 pixels 
per minute in x-direction and by 0.2 pixels per minute in y-direction and then co-added. 



4. Thermophysical modelling on the combined 
dataset 



Our data set is very homogenous with respect to the used 
instrument settings, observing technique, data reduction 
and calibration scheme. The data cover a wide range of 
different observing geometries, including before and after 
opposition, at different wavelengths and rotational phases. 
Therefore, it was possible to adjust the thermal inertia and 
beaming parameters to see how these variations influence 
the calculation of the radiometric diameter and albedo so- 
lutions. The investigations of the resulting diameters and 
albedos give clues about the optimal model parameters. 
Only specific model parameters will allow that e.g. the 
data taken at very different phase angles before and after 
opposition will result in the same radiometric diameter 
and albedo solutions. Another quality criteria is that the 
resulting diameter and albedo values show no trends with 
wavelengths or rotational phase. The goal was to find the 
best diameter/albedo solution with the smallest standard 
deviations which fits all 20 measurements. A similar proce- 
dure was already used by Miiller & Lagerros {HOSI I2002fl 
for several main-belt asteroids. For large, regolith covered 
asteroids, the least-square process gave typical thermal in- 



ertias of 10-15 J m 2 s 0,5 K 1 and beaming parameters of 
p=0.7 and /=0.6 (Miiller et al. lT§55|l . 

As a first step, we tried to find out how the thermal 
inertia influences the determination of the radiometric di- 
ameter/albedo solutions. For observations taken at very 
large phase angles the thermal inertia is the most im- 
portant parameter for a consistent diameter and albedo 
determination. This is due to the non-zero temperature 
of the large terminator which contributes significantly 
to the disk-integrated flux. The thermal inertia T was 
varied in a physically meaningful range between and 
2500 J m~ 2 s~ - 5 K _1 . Where T=0 describes a surface in 
instantaneous equilibrium without any thermal conduc- 
tion into the sub-surface, while r=2500 J m~ 2 s~ - 5 K _1 
corresponds to a highly conductive solid granite sur- 
face without any dust regolith. The Moon, with it thick 
highly insulating dust layer, has a thermal inertia of 
39Jm- 2 s-° 5 K- 1 (KeihmUnH. 

This wide range of thermal inertias has no big im- 
pact on the resulting weighted mean diameter and albedo 
values: A T=0 would give weighted mean values of 
L> cff =0.30km and p v =0.24, while a T=2500 would give 
_D e ff=0.33km and pv=0-19. But the standard deviation 
of the 20 diameter (or albedo) values changes enormously 
(more than a factor of 2) for different thermal inertias. 
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Fig. 2. Equatorial edge-on (top) and pole-on (bottom) im- 
ages of the shape model. 



Figures and 0] illustrate this effect on basis of the 
o 'albedo/ 'albedo and CTdiamctcr/ diameter values. A thermal 
inertia of r=1000 Jm" 2 s~ - 5 K _1 would therefore give 
the best match with our complete observational data set 
(dashed lines). We repeated the whole optimisation pro- 
cess with the highest quality data only (excluding the 
data from April, 9th, 2001, M. Dclbo, priv. comm.). The 
resulting best thermal inertia would then be at around 
750 J m~ 2 s~ - 5 K _1 (solid line). A last robustness check 
with only 9 observations (dotted line) confirmed this so- 
lution. 
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Fig. 3. Thermal inertia optimisation process for the indi- 
vidual TPM albedos and their standard deviation, using 
all 20 individual observations (dashed line), a subset with 
15 observations (solid line) and a subset of 9 observations 
(dotted line). 
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Fig. 4. Thermal inertia optimisation process for the indi- 
vidual TPM diameters and their standard deviation, using 
all 20 individual observations (dashed line), a subset with 
15 observations (solid line) and a subset of 9 observations 
(dotted line). 



In a second iteration, we tested the beaming model, 
parameterised by p, the r.m.s. of the surface slopes, and 
f, the fraction of surface covered by craters. Both parame- 
ters were kept variable between 0.1 and 0.9 (see also Dotto 
et al. I2000[) . However, the effects with p and / are not as 
dramatic, mainly because most of our observations were 
taken at large phase angles where the beaming does not 
play an important role. We could not find a clear minimum 
in the <7 a ibedo/ 'albedo values in the p-f plane. Some good 
solutions disappeared again when we checked for robust- 
ness by using various subsets of the observational data. As 
a conclusion from all optimisation runs we can only say 
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that the very smallest values (0.1-0.3) are very unlikely 
for the beaming parameters p and / for Itokawa, all other 
values still seem to be in agreement with our data set. We 
accepted therefore the default beaming values, p — 0.7 
and / = 0.6, which were derived for main-belt asteroids 
(Miiller et al. 11555) 1 . 

As result of this optimisation process we accepted the 
following values: 

T = 750 J m~ 2 s~ - 5 K _1 , thermal inertia 

p = 0.7 r.m.s. of the surface slopes 

/ = 0.6 fraction of surface covered by craters 

Using the above parameters we derived the weighted 
mean values for D e s and py together with the standard 
deviations 1 : 

D eS = 0.32±0.03km (±0.01 km) 
p Y = 0.19±0.03 (±0.01) 

5. Thermophysical modelling of the new dataset 

We also used the TPM on basis of the established param- 
eters to "transport" or normalise the observed fluxes to 
a reference wavelength of 10.0 pm. Here we concentrated 
on the third dataset which has a 2.5 hour coverage of the 
12.1 h rotation period. For each of the measured values in 
Tbl. we calculated the corresponding D e g and py val- 
ues (for T = 750 J m~ 2 s~ - 5 K~\ ,0=0.7, 7=0.6). This was 
then used again to predict the 10.0 pm brightness for the 
given epoch (see diamond-symbols in Fig. [5)). A weighted 
average D c g and py (of dataset #3) was then taken to pre- 
dict the thermal lightcurve (solid, dashed, dashed-dotted, 
and dotted lines in Fig.EJ), based on the Kaasalainen shape 
and spin- vector model and different thermal inertia values. 

Figure [S] confirms that the thermal inertia has 
to be somewhere in the range between 500 and 
1000 J m~ 2 s~ - 5 K _1 . It also demonstrates that the im- 
plementation of the shape model in combination with the 
spin vector and zero points in time and phase yields con- 
sistent results. Our observational data cover the rotational 
phases between 354.9° (14-Mar-2001 05:50:00 UT) and 
71.1° (Ol-Jul-2004 08:37:00 UT). 

6. Discussion 

All optimisation steps work well under the assumption 
that the shape and spin vector solutions are of good 
quality and that the albedo is the same all over the as- 
teroid surface, e.g. for large main-belt asteroids which 
have almost spherical shapes with a very homogeneous 
albedo distribution due to a thick dust regolith. The shape 
model of Itokawa matches the visual lightcurves taken 
at very large range of phase angles and different observ- 
ing and illumination geometries (Kaasalainen et al. 12005(1 . 
Kaasalainen et al. I(2003|l detected no significant albedo 

1 Accepting the results from the 15 highest quality data only 
decreases the standard deviations significantly. 
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Time [hours], Start time: Ol-Jul-2004 00:03:00 UT 



Fig. 5. Predicted thermal lightcurve at 10.0 pm for the 
time period around the July 1st, 2004 observations. The 
original measurements were "transported" to the 10.0/im 
wavelength via our TPM solution. Predictions and mea- 
surements are shown with their absolute values, no shift- 
ing or scaling in time or flux have been done. 



variegations. Thus, both prerequisites are fulfilled and it 
was for the first time possible to extract thermal prop- 
erties for such an elongated object. So far, this was only 
possible for large main-belt asteroids (e.g. Spencer et al. 
IT555I Miiller & Lagerros lT555)l and a few NEAs (e.g. Harris 
& Davies 11999(1 with almost spherical shapes where the 
radiometric diameter/albedo solutions were not affected 
too much by shape effects or albedo deviations at certain 
epochs. 

Another key element for a successful derivation of ther- 
mal parameters is the wavelength and phase angle cover- 
age of the observational data set. In order to separate 
beaming and thermal inertia effects, it is necessary to 
have equal quality data at different phase angles. An of- 
ten used indicator for reasonable model assumptions is the 
ratio plot "observation/model prediction" (e.g. Miiller & 
Blommaert, r2l70l)l . 

FigureElshows this "obs/mod" -ratio plot for the phase 
angles. Asymmetries in the phase angle plot (see Fig. 
top) confirm the retrograde sense of rotation (see e.g. 
Miiller 2002 for further discussions) and show that the 
terminator has very different temperatures before (a > 0) 
and after opposition (a < 0), resulting in a poor mod- 
elling (large scatter) in the after opposition ratios. The 
TPM requires therefore a higher thermal inertia to match 
the observed fluxes (Fig. bottom) and to eliminate the 
phase angle asymmetry. 

Figure [7| shows how the "obs/mod" -ratio varies with 
wavelengths. In case of a low thermal inertia (top) , one can 
clearly see a trend of the ratio with wavelength. For the 
r=750-case, the scatter of the data points is much smaller 
and the trend with wavelength disappeared. Note, that 
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Fig. 6. The observation/TPM ratios for a thermal inertia 
of (top) and 750 (bottom). The high thermal inertia 
values eliminate the trend with phase angle and reduce 
the scatter significantly. 



Fig. 7. The observation/TPM ratios for a thermal inertia 
of (top) and 750 (bottom). The high thermal inertia 
values eliminate the trend with wavelength and reduce 
the scatter significantly. 



the data points, which are considered to be less reliable 
(Delbo, priv. comm.), are not shown for clarity reasons. 

Our derived albedo of py = O.f 9±0.03 fits nicely within 
the established S-type albedo range of 0.2f±0.07 (Ishiguro 
et al. l2003[) . But this albedo value is strongly connected to 
the Hy-v&lue of 19.9 mag. A change of 0.1 mag in Hy to 
19.8 mag would increase the radiometric albedo by about 
0.02, while the effective diameter would remain practi- 
cally unchanged. Even very large modifications of the Hy- 
value by e.g., 0.5 mag to Hy = 19.4 mag (Hr = 19.0 mag 
from Nishihara et al. 120051 and (V — R) — 0.4 mag from 
Lowry et al. I2005f) would lead to a marginally increased 
D e ff = 0.33 km. But in this case, the directly connected 
albedo would increase to py = 0.29. Based on these con- 
siderations, we give a final solution of py = 0.19^o 03 
to account for the different published Hy values. This 
shows that the size determination is closely coupled to 
the quality of the thermal photometry, while the albedo 
depends much more on the properties from the reflected 
light analysis, i.e., the values in the H-G magnitude sys- 



tem. It should also be noted here that the H-G concept 
is a badly defined convention for irregular bodies. The 
Hy = 19.9 mag is the best possible solution in the context 
of the Kaasalainen et al. (2005) shape model. Nishihara 
et al. I|2005(l also determined a larger G- value of 0.25 (in- 
stead of the 0.21 used here). But this difference would no 
be noticeable in the radiometric results. We also checked 
the robustness of our T-solution against uncertainties in 
the H — G values. But taking the Nishihara et al. (2005) 
H — G values (H = 19.4 mag, G — 0.25) leads to the same, 
very pronounced minimum at around 750 J m~ 2 s~ ' 5 K _1 
in the <7 a ibcdo/ albedo picture (see Fig. OJ). 

The derived radiometric value D c g is the diameter of 
an equal volume sphere, based on the Kaasalainen-shape 
and spin- vector model. A rotating ellipsoid approximation 
can be described by the absolute sizes 2a, 2b and 2c, corre- 
sponding roughly to the x-, y- and z-dimensions of Fig. [5] 
(longest extension, hight in bottom and hight in top im- 
age): 

a = 0.52±0.05km; 2b = 0.27±0.03km; 2c = 0.23±0.02km; 
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Note however, that for very general shapes, like the 
Itokawa-shape by Kaasalainen et al. J2005J), the a/b and 
b/c ratios are not well defined. 

The radar observations (Ostro et al. I2005[l resulted in 
a size estimate of 594x320x288m (±10%), based on the 
same lightcurve-based spin vector by Kaasalainen et al. 
( 2005). The axis ratios of the lightcurve and radar shape 
models agree and have approximate values of a/b=1.9 and 
6/c=l.l. 2 But the effective diameter Z? c ff = 2(afec) 1 / 3 of 
the radar solution is about 15% higher. We tried to use 
the radar effective diameter of about 0.38 km to fit the ob- 
served highest quality data from July 1, 2004. Only with 
very unrealistic assumptions of py—0.5 and a thermal in- 
ertia r=5000 Jm~ 2 s~ a5 K" 1 the TPM predictions would 
provide an acceptable match with the observed fluxes. 
Taking the stated radar uncertainty of ±10% into account, 
we conclude that the radar sizes are overestimated by a 
few percent. 

Sekiguchi et al. pTMjl derived through the NEATM 
(Harris I1998f) radiometric diameter and albedo values 
which agree within the specified error bars with our so- 
lution. The NEATM results for the NEA 2002 NY40 
(data at one phase angle only) also compare well with 
the TPM predictions (Miiller et al. l2tM|l . But for data 
sets covering very different phase angles, as it is the case 
here, the NEATM requires different beaming parameters 
(Delbo et al. 1217113")) . while the TPM can explain all ob- 
served data points with one set of physical and thermal 
parameters. The TPM beaming model, parameterised by 
p and /, can handle the very different illumination ge- 
ometries without artificial correction factors. Additionally, 
the NEATM cannot explain the before/after opposition 
(or morning/evening) effect, unless the beaming parame- 
ters are adjusted differently before and after opposition. A 
detailed comparison with the NEATM was therefore not 
performed. 

We also tried to determine surface roughness proper- 
ties, described in the TPM by p and /. It turned out, that 
the TPM predictions for our given large phase angles are 
almost independent of these beaming parameters. This is 
in agreement with the fact that the effect only plays a role 
at small phase angles where mutual heating within the 
crater structures produce an enhanced amount of infrared 
flux (as compared to a smooth surface). Our phase angle 
coverage therefore does not allow to draw any conclusions 
on the surface roughness, crater structures or r.m.s. values 
of the surface slopes. Additional data at small phase angles 
close to opposition are required for such investigations. 

We also compared the difference between using a 
spherical shape model (together with the true spin vector 
solution under the given observing geometries in combi- 
nation with the TPM) and using the Kaasalainen shape 
model. Assuming a spherical shape gives in fact very simi- 
lar mean (or weighted mean) diameter and albedo values, 

2 Note here again that the dimensions are for a rough and 
not well-defined ellipsoid representation rather than, e.g., for 
the largest extents. 



but the scatter between the 20 derived albedo values is 
about 50% larger. The derived diameter and albedo values 
are then dependent on the rotational phase and to a cer- 
tain extent also on the aspect angle. Additionally, the ef- 
fective diameter and albedo one obtains by simply averag- 
ing our observations is not significantly different than that 
derived from the rotationally resolved observations, since 
they span a range more or less uniformly from maximum 
to minimum of the lightcurve, as shown in Fig. El The 
observations in July 2004 allowed to determine a 10.0 pm 
lightcurve which is perfectly matched by the Kaasalainen 
et al. ( 2005 ) shape and spin vector model together with the 
TPM and a thermal inertia of 750 J m~ 2 s~ - 5 K _1 . In gen- 
eral, measurements taken at relatively large phase angles 
after opposition (here a = —54°), where the terminator 
is still warm (for an object with retrograde rotation), can 
be considered as key observations to determine thermal 
properties of NEAs. The only important point is a suffi- 
cient coverage of the rotational phases. The measurements 
before opposition are much more influenced by the actual 
illumination and observing geometry and only in second 
order by the contribution from the cold terminator. 

The thermal inertia is defined as T — ^/K s p s c s , with 
k s being the thermal conductivity, p s the density and c s 
the heat capacity of the surface material. A dust layer 
on the Itokawa surface, like the one on the Moon, with 
typical Moon-like k s and c s values (Keihm I1984fl would 
require a density several hundred times higher than the 
1250 kg m~ 3 for the Moon to account for the derived ther- 
mal inertia. On the other hand, combining the derived 
high thermal inertia with the S-type bulk density of Ida 
(Belton et al.|l995 ) of 2600 kg m~ 3 and the specific heat 
of Granite c s =800 J kg -1 K _1 the k s value would be in the 
order of 0.3W m _1 K -1 . This seems to be a reasonable 
conductivity for a porous stony material. The porosity 
itself can be determined from the assumed bulk density 
of 2600 kg m~ 3 in combination with an anhydrous ordi- 
nary chondrite surface composition (Ishiguro et al. 2003 ) : 
p = (1 — Pbuik/p) ~ 0.3 (Belton et al. H995fl . In conclusion, 
the mid-IR data are in very good agreement with the as- 
sumption of a bare rock surface. A thick dust regolith can 
be excluded as well as a metallic surface which would have 
a T- value above 10 000 J m~ 2 s~ - 5 K _1 and consequently 
produce a very small thermal lightcurve amplitude (see 
Fig.EJ). 

The total mass of Itokawa is directly connected to the 
assumed bulk density and the determined volume of the 
Kaasalainen-shape model via the D c g value (full uncer- 
tainty range of both quantities has been applied): 
M = pbuik • Volume = 2600 • f 7r(-%a) 3 = 4.5± 2 ;g • 10 10 kg 

7. Conclusion 

The example of Itokawa shows the potential of the TPM 
applications for NEAs. State-of-the-art shape models from 
radar and lightcurve inversion techniques can be used for 
sophisticated thermo-physical investigations. In fact, the 
Itokawa case was the first implementation of Kaasalainen- 
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shape models in the context of the TPM by Lagerros 
lfT^HIITM71HM5|l . The TPM allows the combination of 
observational data taken at different observing and illu- 
mination geometries and wavelengths. No artificial fitting 
parameters are necessary to explain the spectral energy 
distributions nor the thermal behaviour with phase angle. 
And the thermal lightcurve is a " normal" output product 
for objects with known sizes and shapes. 

The Hayabusa mission will characterise Itokawa's 
properties with high reliability. Our derived properties 
can therefore be compared and the TPM be verified. This 
project will establish the "ground truth" for future NEA 
TPM applications. The experience with the thermal data 
are very important for the planning of future observing 
campaigns of NEAs. Depending on the availability of a 
target, the wavelengths, the phase angles and the rota- 
tional phases can be selected in such a way that the ther- 
mophysical characterisation benefits most. E.g. if the sur- 
face roughness properties are of interest, one would have 
to include observations close to opposition. If thermal iner- 
tia is important, certain phase angles and/or a significant 
coverage of the thermal lightcurve are key ingredients for 
a successful study. 
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